Gravitational waves from relativistic rotational core collapse 
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We present results from simulations of axisymmetric relativistic rotational core collapse. The general 
relativistic hydrodynamic equations are formulated in flux-conservative form and solved using a high- 
resolution shock-capturing scheme. The Einstein equations are approximated with a conformally 
flat 3-metric. We use the quadrupole formula to extract waveforms of the gravitational radiation 
emitted during the collapse. A comparison of our results with those of Newtonian simulations shows 
that the wave amplitudes agree within 30%. Surprisingly, in some cases, relativistic effects actually 
diminish the amplitude of the gravitational wave signal. We further find that the parameter range 
of models suffering multiple coherent bounces due to centrifugal forces is considerably smaller than 
in Newtonian simulations. 
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Since more than two decades astrophysicists struggle 
to compute the gravitational wave signal produced by ro- 
tational core collapse supernovae (see e.g. [1]). Besides 
the coalescence of compact black hole and/or neutron 
star binaries, these events are considered as one of the 
most promising (burst) sources to detect gravitational 
waves [2]. Theoretical predictions have been, and still 
are, hampered by two major reasons: Firstly, one does 
not know the rotational state of the core prior to its grav- 
itational collapse, as multi-dimensional evolutionary cal- 
culations of rotating stars have not yet been performed. 
Up to now only quasi-spherical stellar pre-coUapse mod- 
els with a phenomenological description of mixing and 
angular momentum transport are available [3]. Thus, 
predictions from rotational core collapse must rely upon 
extensive parameter studies, i.e., a large set of simula- 
tions. Secondly, reliable collapse simulations require a 
relativistic treatment of the dynamics, as both high den- 
sities and velocities together with strong gravitational 
fields are encountered, as well as the incorporation of 
a realistic equation of state (EoS) and neutrino trans- 
port. All previous studies have either neglected or ap- 
proximated one or all of these requirements. In addition, 
most simulations assumed axisymmetry, the fully three- 
dimensional case still remaining essentially unexplored. 

Concerning the gravitational wave signal the two most 
advanced studies of axisymmetric rotational core collapse 
to date have either investigated a few Newtonian models 
with a realistic EoS and simplified neutrino treatment [4] , 
or a large set of Newtonian models with a simplified EoS 
and no neutrino transport [5]. These simulations pre- 
dict that nonspherical core collapse supernovae emit at 
most ~ IQ-^Mqc^ {Mq = 1.99 x lO^^ g is the sun's 
mass, and c is the speed of light in vacuum) in the form 
of gravitational radiation, with a maximum quadrupole 
wave amplitude A^q ~ 10'^ cm. For a source at a dis- 
tance R this corresponds to a maximum dimensionless 
strain h ~ 10~^^ • (10 Mpc/i?). Similar maximum signal 



strengths have been found in non-axisymmetric Newto- 
nian simulations also using a simplified EoS and neglect- 
ing neutrino transport [6] . 

The present work is a further step to improve the reli- 
ability of the prediction of the gravitational wave signal 
from rotational core collapse. To this end we have per- 
formed the first study (to our knowledge) in the relativis- 
tic framework, extending the axisymmetric Newtonian 
core collapse simulations of Zwerger and Miiller [5]. 

In order to examine a large set of initial conditions we 
approximate the collapsing core by a rotating polytrope 
and use a simplified analytic EoS [7]. We neglect all 
transport effects but include relativistic dynamics in a 
time-dependent spacetime. The pressure P consists of a 
polytropic part Pp and a thermal part Pth, 

P{p, e)=Pp + Pth = Kp^'' + peth(7th - 1), (1) 

where p is the rest-mass density, and £ = £p + eth is 
the specific internal energy which consists of a polytropic 
and a thermal contribution. The polytropic constant K 
and the polytropic index 7p change discontinuously at 
nuclear matter density, p^uc = 2.0 x 10^'* g cm~^, but 
are otherwise assumed to be independent of p elsewhere. 
Pth mimics the thermal pressure in the matter heated up 
by the shock. Writing P as a sum of these two terms 
prevents matter elements from following the same P- 
p-history during infall and expansion after bounce (for 
more details see [7]). In all models we set 7th = 1.5, 
which corresponds to a mixture of relativistic and non- 
relativistic gas. 

Initial configurations are rotating relativistic stellar 
cores in equilibrium computed by means of Hachisu's 
self-consistent field method [8]. All initial models are 
4/3-polytropes, i.e., P = Kp'^^^, with a central density 
of 1.0 X 10^° g cm~^. They rotate according to the differ- 
ential rotation law specified in [8] , in which a parameter 
A (a length scale) determines the angular momentum 
distribution. For large values of A (compared to the size 
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of the initial model) one obtains almost rigidly rotat- 
ing configurations, while small values of A correspond to 
strongly differentially rotating ones. This rotation law is 
the rclativistic extension of the Newtonian one in [9] . The 
second parameter of the initial model is the ratio of the 
amount of rotational energy to the modulus of the grav- 
itational potential energy, /? = ii'rot/l-E'potI- The collapse 
is initiated by suddenly reducing 7p from its initial value 
4/3 to some prescribed value in the range [1.28,1.325]. 
In all simulations we assume both axial and equatorial 
plane symmetry, and do not consider any physical angu- 
lar momentum transport. 

The hydrodynamic evolution of a rclativistic fluid in 
an arbitrary spacetime endowed with a metric g^^, is 
governed by the local conservation laws of the current 
density J'^ = pu^ and the (perfect fluid) stress-energy 
T'"' = phu'^u'^ + Pg^"", where u'^ is the 4- velocity and 
h = 1 + s + P/p is the specific enthalpy. The con- 
servation laws arc the continuity equation V^J'' = 
and the Bianchi identities V ^T'^" = 0, respectively. To 
derive the hydrodynamic equations [10] we introduce a 
set of conserved variables D = pW (rest-mass density). 
Si = phW^Vi (momentum density) and E = phW^ — P 
(total energy density), W = au^ being the rclativis- 
tic Lorentz factor, in terms of the primitive variables 
{p,Vi,s), where Vi is the 3-velocity. The resulting equa- 
tions can be written as a first-order flux-conservative hy- 
perbolic system of conservation laws. 

The time-dependent metric is evolved using the 3-1-1 
splitting of spacetime. In this formalism the Einstein 
equations are written as a system of evolution and con- 
straint equations for the 3-metric 7^ and the extrinsic 
curvature K^^ . Unfortunately, achieving long-term sta- 
ble evolution using the 3 H- 1 system (or improved re- 
formulations) is still an open issue (sec [11] and refer- 
ences therein). Therefore, following [12] we approximate 
the spacetime geometry by assuming that the 3-mctric is 
conformally flat, i.e., = 4>*%j where is a conformal 
factor depending on the coordinates x^, and where jij 
is the fiat 3-nictric. In this approximation, the Einstein 
equations transform into a system of 5 coupled nonlinear 
elliptic equations for the conformal factor (p, the lapse 
function a and the shift vector 

= -2ncl>HphW' - P) - ll^^, (2) 

A{a<P)=2wa<p^ph{3W^~2) + 5P)+'^^^^^^^^^, (3) 

8 

A(3' = mTTa(t>''S' + 2K'^Vj - ^V'Vfe/3^ (4) 

with A and V being the flat space Laplace and Nabla 
operators, respectively. 

The above approximation is exact in spherical symme- 
try. As our focus is to investigate rotational core collapse, 
where the matter distribution does not deviate too much 



from spherical symmetry (contrary to binary neutron star 
mergers [12]), and as the effect of the gravitational radia- 
tion reaction force onto the hydrodynamics is negligible, 
the conformally flat condition (CFC) appears as an ad- 
equate approximation for our purpose. The accuracy of 
the CFC approximation has been considered by Cook 
et al. [13], who find remarkably good results for rapidly 
rotating rclativistic stars in equilibrium, with typical er- 
rors for different variables smaller than 5%. However, the 
quality of the CFC approximation degrades in the case 
of extremely rclativistic nonspherical configurations, e.g., 
rigidly rotating infinitesimally thin disks of dust [14]. 

As off-diagonal elements of •jij are set to zero in the 
CFC approximation, the gravitational wave content of 
the spacetime (the transverse part of the metric tensor) 
is eliminated. Therefore, to compute the gravitational 
wave signal we use the Einstein quadrupole formula in a 
post-processing step. 

In order to exploit the hyperbolic and conservative 
character of the hydrodynamic equations, we have im- 
plemented in our code a Godunov-type flux-conservative 
finite volume scheme based upon the solution of ap- 
proximate Riemann problems at each cell interface (see, 
e.g. [15]). On each time slice, the hydrodynamic data 
are represented by cell averages of the conserved quanti- 
ties D, Si and E. These are propagated to the next time 
level by a conservative Runge-Kutta algorithm, which in- 
volves the numerical fluxes and sources. The metric equa- 
tions (2-4) are discretized on a 2-dimensional grid with a 
9-point stencil. A Newton-Raphson iteration scheme is 
used to obtain the solution for the 5 metric functions. For 
each metric computation, we have to solve a linear sys- 
tem of equations involving the (sparse) Jacobian matrix 
of the discretized equations. 

We have performed a comprehensive number of tests 
of our code. Firstly, we have studied its ability to keep 
rotating neutron stars in equilibrium. To that end we 
have solved the rclativistic hydrostatic equilibrium equa- 
tions for rotating axisymmetric polytropes having differ- 
ent central densities, angular momentum profiles, and 
rotation velocities up to the mass-shedding limit. These 
initial data were evolved both on a static background 
metric and in a fully dynamic spacetime, always find- 
ing stability of the equilibrium configurations over many 
rotation periods. Secondly, we have compared the per- 
formance of our Eulerian code with a May- White-type 
Lagrangian artificial viscosity code [16] in spherical sym- 
metry. When evolving identical initial collapse models 
with both codes we observe agreement in the evolution 
of fluid and metric variables within a few percent. 

In order to pin down the effects of general relativity 
on rotational core collapse, we have compared rclativis- 
tic (R) and Newtonian (N) simulations for a large set of 
initial models. The collapse simulations have been per- 
formed in spherical coordinates using a logarithmic radial 
grid with 200 zones and an equidistant angular grid with 
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TABLE I. Model parameters. Note that model B never 
reaches supranuclear densities. 



Model 


A 


/3r 


/3n 


7p 


7p 




[10» cm] 


[%] 


[%] 


(P < Pnuc) 


(p > Pnuc) 


A 


0.5 


0.47 


0.48 


1.300 


5/3 


B 


1.0 


1.82 


1.81 


1.325 





30 zones. Test runs with more angular zones yield similar 
results. During a typical rclativistic simulation the met- 
ric is calculated every 100th (10th) time step for a central 
density smaller (higher) than O.lpnuc, and extrapolated 
in between. Comparison with a run where the metric 
was calculated every 10th (or every) time step showed 
no significant differences. The total mass / dV^D is 
conserved to better than 10~^ in both relativistic and 
Newtonian rotational collapse. In spherical symmetry, 
the mass is conserved up to machine precision in both 
the Newtonian case, and the relativistic case for a static 
background metric. 

Here we present the results from two representative 
models (Table I). Note that the initial data slightly 
differ with respect to the value of /? between the rela- 
tivistic and Newtonian initial model. The time-evolution 
of the central density is plotted in Fig. 1. Model A 
is a strongly differentially rotating core which collapses 
rapidly (tbounce = 38 ms) due to the large reduction of 7p 
from its initial value 4/3 to 1.30. The model bounces at a 
relatively high density (pc(ibounce) = 9.3 x 10^'' g cm~^) 
because the supranuclear EoS was chosen to be extremely 
soft (7p(p > Pnuc) = 5/3). After bounce the rotating 
core oscillates, in both simulations, with a superposition 
of various radial and axisymmetric modes. As expected, 
the oscillation frequencies of the core arc higher in the 
relativistic model due to the higher densities achieved. 

On the other hand, model B is initially an almost 
rigid rotator, which boimccs later (tbouncc = 98 ms) than 
model A due to its larger subnuclear adiabatic index 
7p = 1.325. The collapse is stopped at subnuclear densi- 
ties (pc(^bouncc) = 3.8 X lO"'^'' g cm"'') by the rapid spin- 
up of the core due to angular momentum conservation. 
After bounce the core first expands, then begins to col- 
lapse once more and rebounccs again. The phenomenon 
of homologous multiple bounces [4] (which occurs with 
a period of roughly 40 ms (50 ms) in the relativistic 
(Newtonian) simulation; see lower panel in Fig. 1) was 
observed in the Newtonian simulations of Zwerger and 
Miiller [5] for strongly differentially rotating configura- 
tions with a sufficiently large initial value of /3 and a 
subnuclear 7p close to 4/3. However, general relativistic 
models show distinct multiple bounces only for extremely 
rapidly and differentially rotating initial models. Most 
of the Newtonian models of Zwerger and Miiller, which 
bounce at subnuclear densities, reach or even exceed in 
the relativistic simulations twice the maximum density 
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FIG. 1. Central density vs. time for model A (upper panel) 
and B (lower panel); relativistic runs are denoted by solid 
lines, Newtonian runs by dotted lines. 



of the Newtonian models. Thus, in many cases the cen- 
tral density largely exceeds nuclear matter density and 
the stiffening of the EoS destroys the coherence of the 
bounce. 

The gravitational wave signature of both models is 
shown in Fig. 2, where we plot the quadrupole wave 
amplitude Afp as defined in [5]. For a source at a dis- 
tance R the corresponding dimensionless strain is = 
8.85 X 10-2^ sin^ 9 ■ (AfoVlOOO) • (10 Mpc/i?), where 6 is 
the angle between the line of sight and the rotational axis 
of the core. 

For the rapidly collapsing model A with its soft 

supranuclear EoS the maximum central density of the 
relativistic simulation is almost 50% higher than in the 
Newtonian case. Nevertheless, the gravitational wave 
amplitude of the relativistic model is lower than that of 
its Newtonian counterpart. This surprising result can be 
understood by noting that in the quadrupole formula the 
gravitational wave amplitude is determined by the sec- 
ond time derivative of the total quadrupole moment of 
the core, which involves a volume integral where the inte- 
grand is proportional to pr**. Thus a higher central den- 
sity due to relativistic effects does not necessarily imply a 
stronger signal. Instead, one must take into account the 
whole density distribution. In Fig. 3 the profiles of p and 
pr'' at the equator [9 = 7r/2) are plotted for both the rel- 
ativistic and the Newtonian model at the time of bounce. 
The density profile of the relativistic model is more com- 
pact than that of the Newtonian model, whereas the den- 
sity of the latter model is larger for radii > 5 km [17]. 
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FIG. 2. Quadrupolc gravitational wave amplitude for 
model A (upper panel) and B (lower panel); relativistic runs 
are denoted by solid lines, Newtonian runs by dotted lines. 



This gives rise to its larger quadrupole moment (as the 
deformation of both models is comparable), and hence to 
its larger gravitational wave amplitude. 

In model B the relativistic effects cause the bounce to 
occur at much higher densities, and the interval between 
the multiple bounces is also shortened. In this case the 
gravitational wave amplitude of the relativistic model is 
actually somewhat larger than in the Newtonian one. As 
we have already pointed out, general relativistic effects 
may play a key role in deciding whether a rapidly rotat- 
ing core bounce will stall at subnuclear density (with a 
very distinctive gravitational waveform as in model B) or 
reach nuclear density (with a waveform as in model A). 
Therefore, relativistic effects have to be taken into ac- 
count when modeling rotational core collapse. We will 
investigate this issue further and present more quantita- 
tive results in an upcoming article. 

As our relativistic simulations yield signal amplitudes 
similar to Newtonian models, the prospects for detection 
of rotational core collapse events are also similar. There- 
fore the predicted signal strength of Galactic events lies 
within the sensitivity range of the upcoming gravitational 
wave interferometers [2]. However, extragalactic events 
seem to be too weak for current detector sensitivities. To 
further improve the waveforms, 3-dimensional effects and 
advanced microphysics have to be taken into account. 

We appreciate contributions and suggestions from 
Y. Eriguchi, J.M- Ibaiiez, and G. Schafer. The calcu- 
lations were carried out at the Rechenzentrum Garching 
and MPI fiir Astrophysik, Garching. 
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FIG. 3. Equatorial profiles of the density p and pr* of 
model A at bounce obtained in a Newtonian (dotted lines) 
and relativistic (solid lines) simulation. 
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